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Abstract 

We study the expected number of real zeros for random linear combinations of 
orthogonal polynomials. It is well known that Kac polynomials, spanned by monomials 
with i.i.d. Gaussian coefficients, have only (2/7r+o(l)) log n expected real zeros in terms 
of the degree n. If the basis is given by the orthonormal polynomials associated with a 
compactly supported Borel measure on the real line, or associated with a Freud weight 
defined on the whole real line, then random linear combinations have n/yj 3 + o(n ) 
expected real zeros. We prove that the same asymptotic relation holds for all random 
orthogonal polynomials on the real line associated with a large class of weights, and give 
local results on the expected number of real zeros. We also show that the counting 
measures of properly scaled zeros of these random polynomials converge weakly to 
either the Ullman distribution or the arcsine distribution. 


1 Introduction 

The expected number of real zeros E[iV n (R)] for random polynomials of the form P n (x) = 
J2k=o Ckxk , wh ere {°k}k =o are independent and identically distributed random variables, 
was studied since the 1930’s. In particular, Bloch and Polya [I] gave an upper bound 
E[7V n (R)] = 0(y/n) for polynomials with coefficients selected from the set {—1,0,1} with 
equal probabilities. Littlewood and Offord I2ZJ-E2B] considered several classes of random 
coefficients, including standard Gaussian and uniformly distributed in {—1,1} or in the 
interval (—1,1). Their results indicated that the expected number of real zeros is actually 
of logarithmic order in terms of n. Shortly thereafter, Kac [20] established the important 
asymptotic result 

E[iV n (R)] = (2/7T + o(l)) logn as n —» oo, 

for polynomials with independent real Gaussian coefficients. In fact, Kac G33-BH found the 
exact formula for E[iV n (M)] in the case of standard real Gaussian coefficients: 


E [N n (R)] 


4 f 1 JA(x)C(x) -B 2 (x) , 

- / -- —- - dx 

7T J 0 A(x) 
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where 


n 


n 


n 


A(x ) = ^x 2 - 7 , B(x) = ^jx 2 - 7 1 and C(x) = ^ j 2 x 2j 2 . 


l=o l=i 1=i 


More precise forms of Kac’s asymptotic for E[Af n (M)] were obtained by many authors, in¬ 
cluding Kac [21], Wang [37], Edelman and Kostlan HU- Wilkins [38] gave an asymptotic 
series expansion for E[7V n (M)] in the case of i.i.d. Gaussian coefficients. 

The asymptotic result for the number of real zeros was generalized by Erdos and Offord 
ra to coefficients with Bernoulli distribution (uniform on {—1,1}), and by Kac [22] to uni¬ 
formly distributed coefficients on [—1,1], Later, Stevens [36] showed that Kac’s asymptotic 
holds for polynomials with coefficients from a certain general class of distributions. Finally, 
Ibragimov and Maslova [18] 19] extended the result to all mean-zero distributions in the 
domain of attraction of the normal law. Many additional references and further directions 
of work on the expected number of real zeros may be found in the books of Bharucha-Reid 
and Sambandham [[[], and of Farahmand [T3] . 

We state a result on the number of real zeros for the random linear combinations of 
rather general functions. It originated in the papers of Kac [29]-[22], who used the monomial 
basis, and was extended to trigonometric polynomials and other bases, see Das i-i and 
Farahmand [13]. We are particularly interested in the bases of orthonormal polynomials, 
which is the case considered by Das [8]. Further generalizations of Kac’s integral formula 
for the expected number of real zeros were obtained by several authors, see e.g. Cramer 
and Leadbetter [7] p. 285]. For any set E C C, we use the notation N n (E ) for the number 
of zeros of random functions fll.ip (or random orthogonal polynomials of degree at most n) 
located in E. The expected number of zeros in E is denoted by E[7V n (i?)], with E[IV n (a, &)] 
being the expected number of zeros in (a, b ) C M. 

Proposition 1.1. Let [a, b] C M, and consider real valued functions gj(x) E (WQa, 6]), j = 
0,... ,n, with go(x) being a nonzero constant. Define the random function 


n 



(l.i) 


1=0 


where the coefficients Cj are i.i.d. random variables with Gaussian distribution Af( 0, a 2 ), a > 
0. If there is M E N such that G' n fx) has at most M zeros in ( a,b ) for all choices of 
coefficients, then the expected number of real zeros of G n (x) in the interval ( a,b ) is given by 



( 1 . 2 ) 


where 


n 


n 


n 


A {x) = ^2g%x), B(x) = ^g j (x)g' j (x) and C(x) = ^[9'^)? ■ (1-3) 


l=o 
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Clearly, the original formula of Kac follows from this proposition for gj(x) = x\ j = 
0,1,..., n. For a sketch of the proof of Proposition 11.11 see [22] • We note that multiple 
zeros are counted only once by the standard convention in all of the above results on real 
zeros. However, the probability of having a multiple zero for a polynomial with Gaussian 
coefficients is equal to 0, so that we have the same result on the expected number of zeros 
regardless whether they are counted with or without multiplicities. 


2 Random orthogonal polynomials 

Let W = e -< ^, where Q : M —> [0, oo) is continuous, and assume that all moments 

/ x j W 2 (x ) dx, j = 0,1, 2,..., 

Jr 


are finite. For n > 0, let 

Pn (x) = Pn ( W 2 , x) = 7 n X n + ... 
denote the nth orthonormal polynomial with y n > 0, so that 

J Pn(x)Pm(x)W 2 (x) dx = j J’ ™ ~ 

Using the orthonormal polynomials {py }y^=o as th e basis, we consider the ensemble of random 
polynomials of the form 

n 

Pn(x ) = y CjPj (x) , n e N, (2.1) 

3=0 

where the coefficients cq, c±, ..., c n are i.i.d. random variables. Such a family is often called 
random orthogonal polynomials. If the coefficients have Gaussian distribution, one can apply 
Proposition [TT]to study the expected number of real zeros of random orthogonal polynomials. 
In particular, Das [SJ considered random Legendre polynomials, and found that E[IV n (—1,1)] 
is asymptotically equal to n/V 3. Wilkins [39] improved the error term in this asymptotic 
relation by showing that E[IV n (—1,1)] = nj\J 3 + o(n £ ) for any £ > 0. For random Jacobi 
polynomials, Das and Bhatt [TO] concluded that E[AT n (—1,1)] is asymptotically equal to nj \/3 
too. They also stated estimates for the expected number of real zeros of random Hermite and 
Laguerre polynomials, but those arguments contain significant gaps. The authors recently 
showed [29] that if the basis is given by the orthonormal polynomials associated to a finite 
Borel measure with compact support on the real line, then random linear combinations have 
n/\/3 + o(n) expected real zeros under mild conditions on the weight. The second and the 
third authors [33] extended this asymptotic to random orthogonal polynomials associated 
with the Freud weights W(x) = e~ c ^ , c > 0, A > 1. The results of this paper provide 
detailed information on the expected number of real zeros for random polynomials associated 
with a large class of weights defined on the whole real line. In particular, they cover the case 
of random Freud polynomials considered in [ 33] . 
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For the orthonormal polynomials {Pjix)}fL 0 , define the reproducing kernel by 


n— 1 

K n (x, y) = ^2p j (x)p j (y), 

3=0 


and the differentiated kernels by 


n—1 

Kn' l) {x,y) = ^2pf\x)Pj\y), 


3=0 


k,l e NU {0}. 


We intend to apply Proposition ll.il with g 3 = pj, so that 


A(x) — K n+ i(x, x), B(x) = Kn+i{x, x) and C(x) — x). (2.2) 

Universality limits for the reproducing kernels of orthogonal polynomials (see Levin and 
Lubinsky |25j-|26j). and asymptotic results on zeros of random polynomials (cf. Pritsker 
[32]) allow us to give asymptotics for the expected number of real zeros for a class of random 
orthogonal polynomials associated with weights from the class ^(C 2 ) (25J. 


Definition 2.1. Let W = , where Q : R —» [0, oo) satisfies the following conditions: 

(a) Q' is continuous in R and Q(0) = 0. 

(b) Q' is non-decreasing in R, and Q" exists in R \ {0}. 

(c) 

lim Q(t) = oo. 

|t| —^OO 


(d) The function 


T{t) 


tQ'{t) 

W’ 


t o, 


is quasi-increasing in (0, oo) 7 in the sense that for some C\ > 0, 


0 < x < y ^ T(x) < C\T(y). 


We assume an analogous restriction for y < x < 0. In addition, we assume that for some 
A > 1, 

Tit) > A in R \ {0}. 

(e) There exists U7 2 > 0 such that 


q '\ x ) < r \< n *)\ 
m*)\ ~ 2 qw ’ 


x eR \ {0}. 


Then we write W e ^(C 2 ). 

Our main result on the global asymptotic for the expected number of real zeros of random 
orthogonal polynomials is below. 
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Theorem 2.2. Let W = G F{C 2 ), where Q is even. If the function T in the definition 
of F(C 2 ) satisfies 

lim T(x) = a G (1, oo], (2.3) 

rr->-oo 

then the expected number of real zeros of random orthogonal polynomials (12. lj) with indepen¬ 
dent real Gaussian coefficients satisfies 


lim -E[7V n (M)] = -L 

n^oo n y/3 


(2.4) 


Theorem 12.21 is a combination of Theorem 12.31 and Corollary 12.51 given below. Define the 
Ullman distribution p, a for 0 < a < oo by 


/40) 


a f 1 t a ~ l , 
K J\x\ Vt 2 - X 2 


X 6 [-1,1], 


and for a = oo, the arcsine distribution by 


Voo( X ) = 


7r 


vT^ 


X- 


X G [—1, 1], 


see LiOj and [24] . We use the contracted version of P n : 


Pn(s) ■■= P n (a n s), n G N, (2.5) 

where a n is the Mhaskar-Rakhmanov-Saff number associated with the weight W, see m, 
[30] . [ 35] and Section 3 below. 

For any set E C C, N*(E) denotes the number of zeros of a random polynomial P*(s) 
located in E. The expected number of zeros of P*(s) in E is given by E[7V*(£!)]. We now 
state the local result on the asymptotic of E[iV* ([a, 5])] for intervals [a, b] C (—1,1). 

Theorem 2.3. Let W = e~® G F(C 2 ), where Q is even. Assume that the function T in the 
definition of F(C 2 ) satisfies (12.31) . If[a,b] C (—1,1) is any closed interval, then 

lim —E [N* ([a, b ])] = ~^=fi a ([a, b ]). (2.6) 

n->oo n yj 3 

We will establish a generalization of Theorem 2.3 for non-even weights in Section 3. Define 
the normalized zero counting measure r n — 1 S Zk for the scaled polynomial Pf (s) of 
(12.51) . where {zfc}^ =1 are its zeros, and S z denotes the unit point mass at z. We determine 
the weak limit of r n for random polynomials with quite general random coefficients {cj}fL 0 . 

Theorem 2.4. Let the coefficients {cj}fL 0 of random orthogonal polynomials (12. ip be com¬ 
plex i.i.d. random variables such that E[| log |c 0 11] < oo. If W = e~ Q G F(C 2 ), where Q is 
even, and if the function T in the definition of F(C 2 ) satisfies ([2.3jl . then the normalized zero 
counting measures r n for the scaled polynomials P*(s) converge weakly to p, a with probability 
one. 
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Related results on the asymptotic zeros distribution of random orthogonal polynomials 
with varying weights were proved by Bloom [5j and Bloom and Levenberg [fTj. Theorem 12.41 
permits us to find asymptotics for the expected number of zeros in various sets. In particular, 
we need the following corollary for the proof of Theorem 12.21 


Corollary 2.5. Suppose that the assumptions of Theorem \2.f\ hold. If E C C is any compact 
set satisfying p a (dE) = 0, then 


lim 1 e K(E) ] = /.„(£). 

n—too n 


(2.7) 


It is of interest to relax conditions on random coefficients Cj, e.g., by considering proba¬ 
bility distributions from the domain of attraction of the normal law as in [T8j EE] . 


3 Proofs 

Our proofs require detailed knowledge of potential theory with external fields, see [21] and 


Let IT be a continuous nonnegative weight function on R such that W is not identically 
zero and lim|, c |_ KX) |x| W(x) = 0. Set Q(x) := — \ogW(x). The weighted equilibrium measure 
pw °f ^ is the unique probability measure with compact support Sw = supp pw C R that 
minimizes the energy functional 


I[u] 


-- dv{t)dv(z) 

x t 


+ 2 


Q du 


amongst all probability measures v with support on R. It satisfies 



dp w (t ) + Q(z) = C, 


Z G S\y, 


and 

[ log t ——- dp w (t) + Q(z) > C, z6l, 

J \z-t\ 

where C is a constant. 

For a weight function W{x) = e~®( x \ where Q is often assumed convex on R, the 
M haskar- Rakhmanov- S aff numb ers 


O—n 0 a n 


are deffiied for n > 1 by the relations 


n = 



xQ'(x) 

\J (*^ n) (fln 


x) 


dx 
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and 


0 = 


Q (z) 


n Ja—n V( X ~ a -n)(a n - x) 


dx. 


We also let 


5 n 2 (®n d~ \ct— n|) aild f3 n c ^( y Q, n d~ n)- 


For even Q , a_ n = — a n , and we may define a n by 


2 

71 


a n tQ'(a n t ) 
) Vi-t 2 


dt = n. 


(3.1) 


Existence and uniqueness of these numbers are established in the monographs IS), (SOI, E3j, 
but go back to earlier work of Mhaskar, Saff, and Rakhmanov. One illustration of their role 
is the Mhaskar-Saff identity: 


which is valid for all polynomials P of degree at most n. We define the Mhaskar-Rakhmanov- 
Saff interval A n as A n := [ a_ n ,a n \. The linear transformation 

. . x — B n _ 

L n {x) = —--, ieK, 

maps A n onto [—1,1]. Its inverse is 

^(s) = (3 n + 8 n S, s G M. 

For e G (0,1), we let 

Jn(^) — d [—1 + e, 1 — e] = [a_ n + e8 n , a n — e8„\. 


The equilibrium density is defined as 


= 


\/(x- a- n )(a n - x) f an Q\s) - Q'(x) 


ds 


71 


s — X 


VTs n) (®n ^ 


, x G A r 


It satisfies the following equilibrium equations [EH p. 41]: 


log -- 7(T n (s) ds + Q(x) = C, x G A n , 


\x — s 


and 


log -j-j-cr n (s) ds + Q(x) > C, x G 


\x — s 


Note that the measure a n (x) dx has total mass n on A n : 

PO'Tl 

/ CT n (x) dx = 71. 
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We also define the normalized version of a n as follows: 


a. 


]0) := W a n( L[ n 1] 0)), s G [-1,1]. 


n 


Note that ct*( s ) ds is a unit measure supported on [—1,1]: 


cFTs) ds = 1. 


'-i 


For details on a n and a*, one should consult the book [2i] . 

In particular, the Ullrnan distribution /i' Q is the normalized equilibrium density for the 
standard Freud weight w(x ) = e~ la ^ a on R, see Theorem 5.1 of [ 35] p. 240], where 

. r(f)T(|) 

7 “ 2r(f + i)’ 

An alternative formula for the Ullrnan distribution follows from that for a n above, namely, 


/ 40 ) = 


2\/l — x 2 f 1 t a — x a dt 

n 2 B a J 0 t 2 - x 2 VI - t 2 


, xe[-l,l], 


where 


B n = 


t° 


7T Jo \J 1 — t 2 
For n > 1, we also define the square root factor 


dt. 


p„(x) = V c - a_„)(a„ - x), x € A, 


(3.2) 


(3.3) 


In the sequel C, Ci, C 2 , • • • denote constants independent of n, x, and polynomials of degree 
< n. The same symbol does not necessarily denote the same constant in different occurrences. 
Given sequences {c n },{d n }, we write 

Cn ~ d n 

if there exist positive constants C\ and C 2 such that for n > 1, 


f 1 1 V c n /d n V (^ 2 * 


Similar notation is used for functions and sequences of functions. 

We start with a general result, our only one that allows non-even weights. I 11 this more 
general setting, P* is given by 


Pn(») = Pn (iUM) • 


rather then by (12.5(1 . 










Theorem 3.1. If W = e Q G IF(C 2 ) and [a,b\ C (—1,1) is any given closed interval, then 
as n —> oo, 

1e re (K 4])] = < +1 (v) dy. 

Proof. The strategy is to apply Theorem 1.6 of (25]. It states that for all r,s > 0, and any 
e G (0,1), we have uniformly for x G J n (e) as n —y oo, 


W 2 (x)K { n’ S \x, X ) 

(cr n (x)) r+s+1 



r-\-s—j—k 

+ °(i)> 


where 

0, j + k odd, 

fc)/ 2 _i_ j j + k even. 

In particular, uniformly in x G J n +i{s), 




WHx)K^fl (x,x) 

cr„+i(x) 

WHx)Rt^(x,x) _ Q'(x) 


( a n+1 (x )) 2 


— 1 + o(l), 

+ °(1)> 


&n+l(x) 


and 


-(i.i) 

n +1 

(cr n+ i(a:)) 3 


W 2 (x)K^{x,x) (Q\x)V 7 T 2 

Un+l^y 3 + ° U ' 


Next, from Proposition 1.1, for any closed interval [l,q\ C J n+ i(e) (where l, q may depend 
on n), 


K n+1 ( x , x ) f x) 


" E[Ar " ([i ’ ?1)1 N ««(*.*) WS ? (*.*>. 

Substituting the above asymptotics, and cancelling, yields 


-E[JV„([/, 9 ])] 

n 


dx. 


-I" 

nir J l 


@n+ 1(x )\ 


n* 


Q'{x) 


1 ) + 


Q'{x) 


, <5i+l (x) J (T rl _)_i(x 

We note that [25} P- 87, Lemma 5.1(a),(d)] uniformly for x G J n+ i(e), 

n + 1 


o(l) + o(l) dx as n —>• oo. 


and 


CTn+l(^) > Cl 

|Q'(x)| <C 2 


Cn+l 

n + 1 


Pn+l(^)’ 
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so that 


Q'(x) 


cr n +i(x) 


< C 3 Sn+1 . < 


C« 


Pn+ 1(^) 7 £(2 — £) 


, X 6 J n+ i(e). 


Thus, uniformly for all intervals [l,q\ C J n+ i(e), as n —> oo, 


-E[/V n ([/,<?])] = — [ a n+ i{x)\l'—+ o(l)dx = 


n 


77.7T 


7T- 


1 + 0(1) n 


n 


73 


o- n +i(x) ete. 


Note that the number N n (E) of real zeros of P n (x) in E equals the number N*(E*) of real 
zeros of P*(s ) in E* := L n (E) = {L n (x) : x G -E}, since L n is a bijection. We recall that a n 
is increasing to +oo and a_ n is decreasing to — oo as n —> oo. It is also known that 

®n+l 1 ,. O—(n+l) .. i ,. ^n+1 

hm -= 1, hm -= 1 and hm —— = 1, 

n->oo a„ n—too a_ iho o ()„ 


see Lemma 3.11(a) of [241 p. 81]. Hence we have 

L n+ i (L^ -1 ^)) = L n+ i(^„ + 5 n s) = ^ 71+1 -+ s as n —>■ oo, 

<Wl <Wl 

uniformly for s in compact subsets of R. If [a, b] C (—1,1), then for large n G N, 

^ ) I®—n + (1 + o ) j ^n (1 7 ] Cl Jn+l (+) j 

provided 0 < e < min{l + a, 1 — b}. It follows that 

1e[at;(M)] = 1e [jv„ (z,J,- 11 ([a,6]))] 

n n 


(3.4) 


l + o(l) 


(n + l)-\/3 JlI 11 (o) 


l + o(l) 
73 

l + o(l) 


■[-i] 


cr n+ i(x) dx 

a* +1 (s) ds (where s = L n+ i(x)) 


(“)) 


r b 

' <r* +1 (s) ds as n —> oo, 


73 Ja ^ 

where we used (13.4p on the last step, and also that 


<+i(7 < 


c 


VT^ 


f, s e ( — I; I); 


□ 


by Theorem l.ll(V) of [24] p. 18]. 

Lemma 3.2. Let W = e~ Q G E(C 2 ), where Q is even. Let a G (1, oo]. If the function T in 
the definition of F{C 2 ) satisfies 

lim T(x) = a, 


then 


lim a n( x ) — f J >'a( X )j x G (—1,1) \ {0}. 
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Remark 3.3. An equivalent form of 


lim T(x) = a G (1, oo) 
£—>•00 


is that uniformly for t in compact subsets of (0,1], 


lim 

£—>•00 


Q'(xt) 

Q'(x) 


= t a ~\ 


Indeed, if this last condition holds, then as x oo, 


T(x)- 1 


Q(x) _ 1 

xQ'(x) xQ'(x) 

Jo Q’{x) 


Q'(u ) du 


t a ~ l dt = 


a 


(3.5) 


Here we also used 0 < Q'(xf)/Q'(x) < 1 and dominated convergence. In the other direction, 
as x —>■ oo, 


Q'(xt ) 
Q'(x ) 


T(xt) Q(xt) T(xt) 
T(x) tQ(x ) tT(x ) 


exp - 


' xt 


Q' 0 ) 

Q(u) 


du 


T(xt ) 
tT(x ) 
T(xt) 
tT(x) 
1 + 0 ( 1 ) 


exp - 


exp 


T(o) 


du 


/xt u 
f x a + o(l) 


du 


' Xt 


U 


exp -(a + o(l)) log - = t a -\ 1 + o(l)). 


Given any e G (0,1), this holds uniformly for t G [e, 1]. 

Proof of Lemma 3.2. We prove the case 1 < a < oo first: 
From (13. ip . as n —* oo, 


n _ 2 r 1 tQ'(a n t) ^ 
a nQ'( a n) ^ Jo Q'(a n )y/ 1 - t 2 
2 7 1 t a 

“t — / / 2 dt = B a . (3.6) 

7T JO yl — t 2 

Indeed, the integrand converges pointwise, and because Q is convex, so Q'(a n t)/Q'(a n ) < 1, 
and we can apply Lebesgue’s Dominated Convergence Theorem. In particular, for n > 1, 
and some C\ > 1 independent of n, 

Cf l n < a n Q\a n ) < Cin. (3.7) 
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Next, we know that for x G (0,1), 


< 0 ) = 


2\/l — x 2 f 1 a n tQ'(a n t ) — a n xQ'(a n x) dt 


n 


n(t 2 — x 2 ) 


vT^I 2 


For t G (0,1) \ {x}, we obtain from (I3.5j) and (13.611 that 


.. ^ntQ J^ntJ 0, n xQ ( d n X ) 
hm - 77 71- 

n->oo n(t z — X z ) 

r a n tQ'(a n t ) - a n xQ'(a n x ) 

= hm - 7=777 w. 2 -7- 

a n Q'{a n )[t z - x z ) 

■j-Gi _ r^OL 

_ P-1 L X 

a f 2 -x 2 ' 


We need a bound on the integrand so as to apply dominated convergence. First, T(it) is 
bounded above. Next, for some £ between t and x, 

&ntQ (ojii) a n xQ (a n x) 
n(t 2 — x 2 ) 

-^(a n uQ (o^/u)) 
n(£ + x) 

< + tfiiQ"( a ntj) 

~ n(t + x) 

Here (13.71) gives (since £/ is increasing) 

&nQ (®n0 ^ &nQ ifln) ^ Cl 

n(f + x) — n(t + x) ~ x 

By definition of ^(C 2 ) and boundedness of T, we have 


so that 


n(t + x) 

Thus, for all t G (0,1), 


Q"(y) CaT (y) C 3 

0 , J/>0, 

Q(y) v y 


a n €Q"(an€) / ^ o 2 ^Q'(«nO / ^ a n Q'(a n £) C 4 

— L '3- y. , - 7 S G 3 -;-— -. 


a n £n(t + x) 
u n ((V (^62^) a n xQ ( a n x ) 


n(t + x) x 


n(t 2 — x 2 ) 


< 


CV 


x 


and we can apply dominated convergence to deduce that 


lim <(x) 

n—>• 00 


2\/l — x 2 f 1 t a — x“ dt , 

Jo = 


12 

































Next, we deal with the case a = oo: 

Let 0 < r < s < 1. We consider x G (0, r] and split 


2VT 


x* 


71 

= :/i + / 2 - 


+ 


a n tQ'(a n t ) — a n xQ'(a n x) dt 


Tl(t 2 — x 2 ) 


Vl -t 2 


(3.8) 


We shall show that the main contribution to cr* comes from 1-2■ Since the integrand in the 
integral defining <7* is nonnegative, we have for x G (0, r] that 


U = 


< 


< 


< 


2y/l ~ X 2 

r 

TT 2 J 

L 

2y/l - X 2 

f 1 

n 2 J 

i 

s 

2y/l^X 2 


7r 2 (s 2 — x 2 ) 

\n t 

y/1 — X 2 

2 

7r(s 2 — x 2 )n 7r 

y/1 — X 2 



n(t 2 — x 2 ) 


vT^t 2 


n{t 2 - X 2 ) y/1 - t 2 

f 1 \ dt 

/ a n tQ ( a n t ) 


OjntQ 


vT^t 2 

dt 




7r(s 2 — X 2 ) 

Next, note that by the lower bound in (3.5) of [23] p. 64], for t G [0,r], 

Q < antQ'(a n t ) < a n rQ'(a n r ) < T(a n r ) fr^™™{^c 6 T(a n r)} 


Oj n sQ (u n s) a n sQ ( CL n s ) T(o, n s^ \s 

/ r\ CsT(a n r) 

^(;) ' 

since T is quasi-increasing. Our hypothesis 

lim T(x) = oo 

x—>oo 

gives 

,. dntQ (Qn^) n 

hm max —— = (J. 

n—>oo te[o,r] a n sQ'(a n s) 

It also then follows easily from (13. ip that for each fixed r G (0,1), 

,. a n rQ'{a n T) 

hm -= 0. 

n—>oo ft 


(3.9) 


(3.10) 


(3.11) 
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Now uniformly for x G [0, r], 


h > 


> 


n 


2\J\ — x 2 f 1 a n tQ'(a n t ) — a n xQ'(a n x) dt 

vr 2 (l -x 2 ) J s 

1 2 r 1 

7T\/1 - X 2 7m J s 
l + o(l) 2 f 1 


a n tQ'(a n t)( 1 + o(l)) 


x/T^t 2 

dt 




TTy/l - X 2 7m ,/ 0 
1 + 0 ( 1 ) 


a n tQ {flrd) 


dt 


Vi-t 2 


7T\/1 — X 2 


as n —> oo, 


(3.12) 


by (13. ip and using (13.101) . Now we deal with A - it clearly suffices to show only an upper 
bound. Let s < p < 1. By definition of the class FiC 2 ) and (13.lip . we have that 


h = 


< 


2a/1 — x 2 f s a n tQ'(a n t ) — a n xQ\a n x) dt 


7T 


2Vl zr x 2 


- max 

7 T z nX uG[0,s] 


d_ 

du 


n(t 2 — x 2 ) 
a n uQ (yOj n xi^ 


VT^ 

" s dt 


1 0 


vT^I 2 


Co 


< — \a n Q'(a n s) + max atuQ"(a n u)] 

TlX wG[0,s] 

C 

< o(l) H-- max a n Q'(a n u)T(a n u) as n —> oo. 

TlX uG[0,s] 

Using the fact that T is quasi-increasing and the lower bound in (3.5) of [23, p. 64], we 
continue this as 


C 

h < o(l) H- -a n Q'{a n s)T{a n s) 

nx 

^ / -. \ | Cg ,, /T(a n s) fs 

— ^( 1 ) UnQ V^nPjTjT/ 7 ( — 
nx T[a n p) \p 

C (s 

< o(l) H- -a n Q'{a n p) sup ( - 

2/G[0,oo) \P 


max{A,C6T(a n s)}—1 


T(a n s ) 


max{A,C6?/}—1 


y = o(l) as n —* oo, 


by (13. lip and as s/p < 1. Together with the fact that I\ > 0, and using 
we have shown that for x G (0, r], 


, dam), 


< lim inf a* (x) < lim sup a* (x) < 


VT 


X z 


TTy/1 — X 2 n ^°° ny ' n.-> oo 1 " nV 7 7r(s 2 37 2 ) 

As s is independent of r, we can let s —> 1— to deduce that for x G (0, r], 

1 , , . 


lim a*(x) = ■ 
n—>oo 7T\/1 — X* 


= PooK X ) 


□ 
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Proof of Theorem II?. 31 We know from Theorem 13.11 that 


iE[J\C([o,6])] = < + Mdy- 

Lemma [3.21 gives for 1 < a < oo that 

hm cr* n+1 (y) = fi a (y), y G (-1,1) \ {0}. 

n—>oo 

Next, by Theorem l.ll(V) of [2H p. 18], 


<+i0) < 


C 


VT^ 


p s G (-1,1). 


Lebesgue’s Dominated Convergence Theorem now implies that 


linr -E [N* ([a, b})} = — / hm a* n+1 (y) dy = -r=/x a ([a, b ]). 
n-xx>n V3 Ja n ^°° V3 


Lemma 3.4. IfW = e ® G F(C 2 ) then 


hm al /n = 1. 

n—>oo 


□ 


Proof. Lemma 3.5(c) of [24] p. 72] implies that there is a constant C > 0 such that 

1 < — < Cn 1 ^ for all n > 1, 

a i 

which immediately gives the needed result. □ 

Lemma 3.5. Let W = e~ Q G F(C 2 ), where Q is even. If the coefficients {cj}°L 0 of random 
orthogonal polynomials (12.ip are complex i.i.d. random variables such that E[| log |co||] < oo, 
then 

hm \\PnWW 1 - = 1 with probability one. 

Proof. Using orthogonality, we obtain for polynomials defined in (12. ip that 


Hence 


P n (x)\ 2 W 2 (x ) dx = ^ |Cj| 2 

3=0 


( [•oo \ l / 2 

/ \P n (x)\ 2 W 2 (x) dx I < (n + 1) max |c,- 

d-oo ) o<j<n 
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Lemma 4.2 of 


lim 

n—¥ oo 


see (4.6) there) implies that 

, V(2n) 

|P n (x)| 2 hh 2 (a;) dx 


= lim ( 


max \Cn 


\o<j< 


n 


1/n 


= l 


with probability one. That is, 

lim llPnVLlI 1 /” = 1 with probability one. 

n—>-oo ^ W 

We use the Nikolskii inequalities of Theorem 10.3 of [24J p. 295] stated as 


(3.13) 


law'll, 


< C, 


n 


1/2 


(r(a„)) 1 / 4 ||P„W'|| I1(i , | 


and 


l|f„»'|| 11(E) < c 2 a ]! 2 . 

Since T(a n ) = 0{n 2 ) by Lemma 3.7 of [24] p. 76], we obtain that 

±-^=\\P n W\\ 

L 2 (M) — \\PnW\\ Lao{R) <C 3 n\\P n W\\ L2{R) , 

^2 V Un. 


□ 


and the result follows by applying Lemma [3741 and (13.131) . 

Lemma 3.6. Let W = G IF(C 2 ), where Q is even. If the function T in the definition 
of IF(C 2 ) satisfies 

lim T(x) = oo, 


then 


lim 7 l /n a n = 2 , 


where is the leading coefficient of the orthonormal polynomial p n (x) associated with the 
weight W 2 . 

Proof. Theorem 1.22 of [24] p. 25] gives 


1 f 

7n = \2 


-n-l i r an , Q(s) ds 

10 2 7r J —cin 

e 


—. v y ds 

V4 -« 2 (i + o(l)) as n —> oo, 


so that 

1 j-a-n Q(s) 

V I T17T J — a n / 7 17 

:ii n a n = 2a n 2n e V“™- s (1 + o(l)) as n —> oo. 

Since Q is increasing on (0, oo), we have that 


1 f a 

0 < lim — 

n-K>o riTT 


-a n \J { 


Q(s) Q(a n ) f cln ds Q(a n ) 

ds < - — 


a l - s 2 


nir 


I-On V a n - s2 


< 


c 


n 


y/T(an) 


(3.14) 


->■ 0 (3.15) 


as n —> oo, by Lemma 3.4 of [24] p. 69]. Thus 


1 

lim — 

n—>oo TITT 


Q(s) 


—a„ \/i 


a 2 n - s 2 


ds = 0 , 


and (13.14ft together with Lemma [3.41 imply the result. 


□ 
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Lemma 3.7. Let W = e~® G IF(C 2 ), where Q is even. If the function T in the definition 
of IF(C 2 ) satisfies 

lim T(x) = a G (1, oo), 

£—>•00 

then 

lim 7 \l n a n = 2 e 1/o , 

n— 

where 7 n is the leading coefficient of the orthonormal polynomial p n (x) associated with the 
weight W 2 . 


Proof. Considering Lemma [3.41 and (I3.14p . we only need to show 

lim — r Q(s) ds = lim I [' Q M dt = If 

n ~+00 nn J~ an ^a? n - s 2 n ^°° n J_ 1 7r\/l - t 2 

In terms of the function T, we can recast this as 

.. If 1 1 a n tQ'(a n t ) 

hm — / —--- . dt = la. 

n^oo n J_ x T{a n t ) 7 T \/1 — t 2 


Using our assumption that lim^ 00 T(t) = a G (l,oo), we have uniformly for \t\ > a n 1 ^ 2 , 
that T(a n t ) = a(l + o(l)), so as the integrand is non-negative, 


a n tQ'(a n t ) ^ _ 1 + o(l) 1 


a n tQ (^n^) 

n Ja~ 1/2 <\t\<l T(a n t ) 7T\/1 - t 2 « n 7a“ 1/2 <|i|<l 7T\/l - t 2 

The integral over the remaining range is small: for j = 0,1, using 

1 a n tQ (a n t) 


dt. 


(3.16) 


and Hindoo a„ = 00 , 


1 

0 < - 


n J\t\<a~ 1/2 T(a n t)i nVl - t 2 


dt 


Tl< 

< 1- al/ 2 Q'(al/ 2 ) 


n 


1 Ahr y/1 - a n 1 
Thus (13.161) and (13.ip yield 


2 < 1/2 < C 


QPP <C Q K) (1) 


n 


n 


antQ (a n t) 


n T(a n t) TTy/l - t 2 


dt = 


1 + o(l) 1 f l a n tQ'(a n t) 1 + o(l) 


a 


n J 7r\/l — t 2 


dt = 


a 


□ 


Proof of Theorem 2.4 


We first deal with the case 


hm T(a;) = 00 , 
£—>•00 
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and show that the normalized zero counting measures r n for the scaled polynomials P*(s) 
converge weakly to the arcsine distribution with probability one. Theorem 2.1 of 0 p. 
310] states that if {M n }“is any sequence of monic polynomials of degree deg(M„) = n 
satisfying 

lim sup (3.17) 

n—>oo ^ 

then the normalized zero counting measures r n for the polynomials M n converge weakly to 
/loo. Note that 1/2 in the above equation is the logarithmic capacity of [—1,1], see Corollary 
5.2.4 of [34] p. 134]. We show that the monic polynomials 


M n (x) \= P*(x)/(c n ^ n a^), n6N, 

satisfy (13 .1 71) with probability one, so that the result of Theorem 12.41 follows for a = oo. We 
know from Lemma [3.51 that 

lim sup < 1 with probability one. 

n— >og 


Using the contracted weight 

w n (s) \= y/W(a n s ) = e s G M, 

and the properties of a n m p. 4], we obtain that 

II II L°°([—1,1]) = || II L°°([—a n ,a n ]) = II PftW || L°°(R) ' 


It follows that 


lim sup ||P/w”||]//](r_ 1 n) < 1 with probability one. 

n—>oo ’ 

Since lim n ^ooQ(a n )/n = 0 (recall (13. 1 5jl ) . we have that 

lim snn II P*|| 1//n < lim snn II P* 7/; n ll^ n P Q(on)/n < i 

mil sup ||^ n || i oo(r _ 1 m ^ 11111 sup || r n UJ n\\L co ([-l.l\) e — 1 

n—>• oo n—>-oo 


with probability one. We use below that lim^,*, ^(n n a n = 2 by Lemma 13.61 and that 
lim^oo | c n | 1,/n = 1 with probability one by Lemma 4.2 of [[32]. This implies that 


lim sup ||M n ||j£ ([ _ 1;1]) 

n—>• oo 


lim sup 

n—>-oo 


p* 


1/n 


1/n 

L°°([-l,l]) 

1 


lim sup ||P nl lL o ° ([ - 1 ,i ]) TU TT/r7 1/7 

n—ioo \ c n\ ' y„ 



with probability one. 

Next, we prove the case 


lim T(x) = a G (1, oo). 

x —>-oo 
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Recall that the standard Freud weight with index a is given by 

w(s) = e _7 “ |s| “, seR, 

where 

r(f)r(|) f 1 t°-' 

7 “ 2 r(f + i) X Vtt 

see [35| p. 239]. Since ") a +i = B a n/2, we apply r(l/2) = s/tt and T(t + 1) = tT(t) to obtain 
that 

„ _ 27 ^ 2 r(f)r(i) r(s±i)r(|) 1 

7 “ “ 7 “ IT x2r(f + i)2r(3±l + i) a 

Note that by [321 P- 240], F w = log 2 + 1/a is the modified Robin constant and n w = 
is the equilibrium measure corresponding to w. Following [35], we call a sequence of monic 
polynomials {M n }™ =1 , with deg(M n ) = n, asymptotically extremal with respect to the weight 
w if it satisfies 

lim JKMJI/" = e~ Fw = e~ l ^ a /2. (3.18) 

Theorem 4.2 of [35] p. 170] states that asymptotically extremal monic polynomials have 
their zeros distributed according to the measure fi w . Namely, the normalized zero counting 
measures of M n converge weakly to u w = u a . On the other hand, by Corollary 2.6 of [321 p. 
157] and Theorem 5.1 of [321 P- 240], 

|| ||ioo([—ip]). 

Together with Theorem 3.6 of [35] p. 46], (13.181) is equivalent to 

lim sup < e~ Fw = e~ l,a /2. 

n—too ’ 

We show that the monic polynomials 

M n (x) := P*(x)/(cn 7n<), n G N, 

are asymptotically extremal in this sense with probability one, so that the result of Theorem 
12.41 follows. Note that 

lim WPM)'- = 1 with probability one 

n—>-oo ^ W 

by Lemma 13.51 and that 

ll^n W rJlz,°°([—1,1]) = ll^^ / HL°°([-a r[ ,a r[ ]) = 11 W \ \ L oo ( R ) 

by [2H p. 4], Hence 

lim sup ||P :<|!^ ([ _ 1 < 1 with probability one. 

n—>oo 
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By Lemma [3.71 and since linage*, |c n | 1//?l = 1 with probability one by Lemma 4.2 of (32] . it 
follows that 


limsup \\M n w n \\ 1 /™n 1>1]} 

n—too 


= limsnp ||P> n ||i / » ([ _ 1)1]) 1/n \ /n 

Cn 7 n a 

= ^ limSU Pll P > n |lL“([-l,l]) 

n—>-oo 

= e~ Fw limsup ||i 3 *ty”||i / » f[ _ 1 in • 


n 


On the other hand, 

limsup \\P*w n \\ l lZn_ uj) < limsup ||P,>”|| P ” ([ _ 1A]) \\w/w n \\ Loo{[ _ m 

n—>-oo ’ n—> oo 

< limsup ||w/wn|| L oo([_i ) i]) ■ 

n—> oo 

Since and w are both even, it remains to show that 


limsup \\w/w n \\ L oo ([0jl]) < 1- 


Let £ G (0,1). For x G [e, 1], (13. 6 p and then (j3.5[) give that 
Q(a n s) _ 1 + o(l) f s a n Q'(a n x) _ 1 + o(l) 


n 


B n 


1 0 &nQ (®n) 


B, 


x a -\l + o{l))dx 


a JO 


aB, 


-(1 + o(l)) = 7 Q s a (l + o(l)) as n —* oo. 


This holds uniformly for s G [e, 1] as (13. 5 p does. Hence 

Q (Ora^) 


I / II ( ( 

\ W / W n\\L°°([£,l]) = SU P eX P ‘ 
s6[e,l] V 


n 


— 7 a s —> 1 as n —y oo. 


Since Q is increasing, we also have that 


| w/w. 


T*-llz,°°([ 0 ,e]) 


< exp 


Q(a n e ) 


n 


-> exp(7 Q £“). 


We finish the proof by letting e —y 0. □ 

Proof of Corollary QQ1 Consider the normalized zero counting measure r n = - Y^k =i <^z fc for 
the scaled polynomial P*(s) of (12.51) . where {zk}f. =l are the zeros of that polynomial, and 
S z denotes the unit point mass at z. Theorem 12.41 implies that the measures r n converge 
weakly to fi a with probability one. Since / i a {dE ) = 0, we obtain that t u \e converges weakly 
to /j, q \e with probability one by Theorem 0.5' of (23] and Theorem 2.1 of |2J. In particular, 
we have that the random variables r n {E) converge to /i a (E) with probability one. Hence 
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this convergence holds in L p sense by the Dominated Convergence Theorem, as r n (E) are 
uniformly bounded by 1, see Chapter 5 of HU. It follows that 

lirn E[| r n (E) - n a (E) |] = 0 

n—>oo 

for any compact set E such that /i a {dE) = 0, and 

\E[T n (E) - Ha{E)] I < E[| r n (E) - Ha(E) |] —> 0 as n —>■ oo. 

But E [r n (E)\ = K[N*(E)]/n and E [p a (E)\ = /j, a (E), which immediately gives (12.711 . □ 

Proof of Theorem 12.21 Theorem 12.31 gives that 

lirn -E [N* ([a, 6])] = -^ Q ([a, b}) 

n->oo n v3 

for any interval [a, b\ C (—1,1). Note that both E [IV* ( H )] and are additive functions 

of the set H. Moreover, they both vanish when H is a single point by (|2.7j) and the absolute 
continuity of with respect to Lebesgue measure on [—1,1], Hence (12.7p gives that 

lim i-E [N‘ (R \ (-1,1))] = MR \ (-1.1)) = 0. 

n— >oo n 

It now follows that 

lim —E[JV^(R)] = -h^,((-l, 1)) = -h. 

n^oo n -y/3 v 3 

To complete the proof, observe that IV*(M) = IV n (M), so that E[IV*(M)] = E[IV n (M)], since 
L n (x) = x/a n is a bijection for each fixed n. Therefore (12.41) is proved. □ 
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